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This article investigates the role of the Saharan Air Layer (SAL) in tropical cy- 
clogenetic processes associated with a non-developing and a developing African 
easterly wave observed during the Special Observation Period (SOP-3) phase of 
the 2006 NASA African Monsoon Multidisciplinary Analyses (NAMMA). The 
two waves are chosen because both interact heavily with Saharan air. A global 
data assimilation and forecast system, the NASA GEOS-5, is being run to pro- 
duce a set of high-quality global analyses, inclusive of all observations used op- 
erationally but with denser satellite information. In particular, following previous 
works by the same Authors, the quality-controlled data from the Atmospheric In- 
frared Sounder (AIRS) used to produce these analyses have a better coverage than 
the one adopted by operational centers. From these improved analyses, two sets of 
3 1 5-day high resolution forecasts, at horizontal resolutions of both half and quar- 
ter degrees, are produced. Results show that very steep moisture gradients are 
associated with the SAL in forecasts and analyses even at great distance from the 
Sahara. In addition, a thermal dipole (warm above, cool below) is present in the 
non-developing case. Moderate Resolution Imaging Spectroradiometer (MODIS) 
show that aerosol optical thickness is higher in the non-developing case. Alto- 
gether, results suggest that radiative effect of dust may play some role in produc- 
ing a thermal structure less favorable to cyclogenesis. Results also indicate that 
only global horizontal resolutions on the order of 20-30 kilometers can capture 
the large-scale transport and the fine thermal structure of the SAL, inclusive of the 
sharp moisture gradients, reproducing the effect of tropical cyclone suppression 
which has been hypothesized by previous authors from observational and regional 



modeling perspectives. These effects cannot be fully represented at lower resolu- 
tions. Global resolution of a quarter of a degree is a minimum critical threshold 
to investigate Atlantic tropical cyclogenesis from a global modeling perspective. 



1. Introduction 


African Easterly Waves (AEWs) have been recognized as prominent weather- 
producing events of northern tropical Africa (e.g Burpee 1974; Asnani 2005) and 
have been extensively studied from observational and modeling perspectives (e.g 
Hsieh and Cook 2005; Kiladis et al. 2006). However, the development of AEWs 
into tropical depressions remains one of the most challenging problems in the 
prediction and modeling of Atlantic tropical cyclones. 

Investigations of tropical cyclogenesis have been generally carried out from 
either strictly observational or high-resolution mesoscale points of view, because 
the lower resolution used in most global models is deemed to be inadequate to 
trigger spontaneous cyclogenesis. However, in recent years a number of global 
models have reached the resolution of 10-40 km and have started to display some 
tropical cyclogenesis capability. Atlas et al. (2005) shown that a semi-operational 
version of the model being run at NASA at that time was capable of produc- 
ing well-defined tropical cyclone structures even if initialized with low-resolution 
analyses in which cyclones were either not present or analyzed properly. Shen et 
al. (2006) used a higher resolution research version of the same model producing 
center pressures on the order of 920 hPa for hurricane Katrina; Reale et al. (2007) 
analyzed the performance of the European Center for Medium Range Weather 
Forecasts (ECMWF) so-called Nature Run (a 13 -month long simulation in free 
running mode performed with the then-operational ECMWF model) and verified 
that the simulation contained realistic AEW activity and a realistic number of At- 
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lantic tropical cyclones. 

Strictly speaking, the resolutions adopted by the above referenced studies can- 
not be considered adequate yet to resolve the fine structure of tropical cyclogenetic 
processes, because the organization of several convective centers into a rotating 
system requires cloud-resolving resolutions which are not yet possible in real- 
time global numerical weather prediction models. Despite these limitations, it is 
nonetheless remarkable that the aforementioned high-resolution global models, 
with different experimental or operational configurations, display distinct cyclo- 
genetic capabilities. 

Moreover, it is important to emphasize that the processes of cyclogenesis and 
cyclogenesis suppression should be studied from different perspectives, and global 
models do have the unquestionable advantages of 1) better capturing the large- 
scale forcings involved, and 2) not relying upon somewhat questionable boundary 
conditions imposed on the domain’s boundaries of limited area models. 

Among the large-scale problems connected with tropical cyclogenesis, one of 
the most debated is perhaps the role of the Saharan Air Layer (SAL). The SAL, 
a layer of hot dry air rich in dust and produced over the Sahara Desert (Carlson 
and Prospero 1972), has been subject to intense scrutiny in order to understand its 
possible connection with tropical cyclone development or suppression. 

Dunion and Velden (2004) investigate the SAL with the aid of Geostationary 
Operational Environmental Satellite (GOES) and recognize it as a possible mech- 
anism linked to tropical cyclone suppression, emphasizing three possible cycloge- 
nesis inhibiting actions: increased static stability, increased vertical shear and dry 
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air inclusion. On a longer temporal scale, Lau and Kim (2007) present evidence 
on that the mean activity of the SAL is likely to impact the large-scale environ- 
ment, including sea-surface temperature (SST) and mean vertical shear, during 
the entire hurricane season. 

However, Jenkins and Pratt (2008) have investigated the connection between 
dust, observed lightning and tropical cyclones during NAMMA-06, with empha- 
sis on mycrophysical processes, confirming that SAL increases static stability and 
may reduce SST, but also suggesting that SAL may invigorate rain bands con- 
nected with cyclogenesis, and can also increase cyclonic vorticity of AEWs. In 
the subsequent paper by the same first Author, Jenkins et al. (2008) present more 
evidence of a possible tropical cyclone rain band invigoration due to Saharan dust. 

At this time, it appears that there is no conclusive evidence on a clear, unam- 
biguous role of the SAL with respect to cyclogenesis, and the only certain fact is 
the extreme complexity of the interaction between the SAL and tropical weather 
systems, which requires much more research. 

In this paper, it is shown that a high resolution global data assimilation system 
(DAS) coupled with a global model of comparable or even higher resolution are 
suitable tools to investigate the role of the SAL not only on the wave scale, but 
also from the point of view of the large-scale transport from its source region, and 
can therefore represent the possible modifications of SAL as the waves propagate 
over thousands of kilometers. 

At the same time, it is shown that model horizontal resolution not lower than 
a quarter of a degree is needed to unveil some of the SAL thermodynamic and 
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kinematic features, such as the increasingly narrow structure of the dry air fila- 
ments being intruded in a tropical circulation, and the sharpness of the boundaries 
between Saharan and non- Saharan air. 

In this work, we explore processes associated with two waves which were ob- 
served during the Special Observing Phase (SOP-3) of NAMMA, with the aid 
of the global data assimilation and forecasting system called NASA GEOS-5. Of 
particular concern for this article is the representation of the thermal structure dur- 
ing cyclogenesis, and in particular the possible role of SAL. Section 2 describes 
the data assimilation, computational aspects and simulations performed, Section 
3 analyzes the results, Section 4 presents a general discussion and Section 5 states 
the conclusion of this work. 

2. The NASA GEOS-5 Data Assimilation and Fore- 
casting System 

The effect of Saharan air on Atlantic tropical cyclogenesis is investigated using 
a high-resolution global DAS based upon the Gridpoint Statistical Interpolation 
(GSI) analysis algorithm co-developed by the National Centers for Environmen- 
tal Predictions (NCEP) Environmental Modeling Center, documented in Wu et al. 
(2002). The NASA GEOS-5 combines the above referenced DAS with the NASA 
atmospheric global forecast model documented in Bosilovich et al. (2006) and 
used, among others, by Reale et al. (2008) and Reale et al. (2009). The fore- 
cast model shares the same dynamical core (Lin, 2004) with the so-called NASA 
finite-volume General Circulation Model (fvGCM) which, in the version known 
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as GEOS-4, demonstrated remarkable capabilities in hurricane forecasting (Atlas 
et al., 2005; Shen et al., 2006). 

The GEOS-5 however contains a new physics (convective and boundary layer 
parametrizations), developed predominantly by the Global Modeling and Assim- 
ilation Office (GMAO). The changes in physics made the system quite different 
from its predecessor. From one side global forecasting skill has substantially in- 
creased, with the performance of the model becoming comparable to major oper- 
ational centers (e.g. Reale et al. 2008). Moreover, the heavy reliance on satellite 
data has made the performance quite symmetric between the two hemispheres, un- 
like many operational models which rely more on conventional observations and 
consequently display a substantial difference in performance between the north- 
ern and the southern Hemisphere. However, the capability of representing very 
deep tropical cyclones, which was remarkable in the previous versions (Atlas et 
al. 2005; Shen et al. 2006) seem to have disappeared in the current configuration, 
possibly because of a different convection parametrization adopted. 

However, the most important difference between the GEOS-5 and the sys- 
tem used in these two previous studies is in the initialization. Both Atlas et al. 
(2005) and Shen et al. (2006), despite running the model at horizontal resolu- 
tions of 0.25° and 0.125° respectively, used initial conditions at about 1° resolu- 
tion obtained from the NCEP operational analyses. On the other contrary, this 
study makes full use of the data assimilation system included in the GEOS-5, and 
produces its own set of global, high-resolution 0.5°, high-quality analyses. The 
analyses contain all the operational data used by NCEP at that time, but a much 
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denser satellite information. In particular, a much improved coverage from the 
Atmospheric Infrared Sounder (AIRS) is utilized, as compared to the one used 
by operational forecasts. Instead of clear-sky AIRS radiances, which make use of 
AIRS data only from channels completely unaffected by clouds, we use quality- 
controlled temperature retrievals obtained also under partly cloudy conditions, 
following Susskind et al. (2006) and Susskind (2007). Despite the widespread 
general assumption that clear-sky radiances are the best way to assimilate AIRS 
data, Reale et al (2009) have shown retrievals under partial cloudy conditions to 
provide better analyses of a tropical cyclone in the Indian ocean than the clear-sky 
radiance methodology. Finally it should be stressed that the DAS for this exper- 
iment is configured without any bogus vortex or vortex relocation technique. As 
a consequence, tropical development, if any, is spontaneously produced by the 
model. 

One 35-day global data assimilation, starting at OOz 13 August 2006, is per- 
formed with the GEOS -5 DAS, to cover the entire period of the SOP-3 phase of 
the NAMMA campaign. From the analyses, two sets of 31 five-day forecasts, at 
two different resolutions, are performed. The DAS and forecast model configura- 
tion is the same used by Reale et al. (2008) but at higher resolution: the DAS is 
being run at a horizontal resolution of 0.67° x 0.5° and the forecast model is run 
at and 0.25° x 0.33° and 0.67° x 0.5°, both with 72 vertical levels. 

Analyses and forecasts are validated and verified against operational NCEP 
analyses. The analyses correspond remarkably well to the operational NCEP 
analyses, as it will shown later, but with more information available, due to higher 
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resolution and possibly denser satellite information assimilated. Global and hemi- 
spheric forecast skills, investigated through the anomaly correlation for 500 hPa 
geopotential as a function of forecast time, range from approximately 0.76 to 0.80 
at day 5 (not shown), indicating that the forecast skill for boreal summer condi- 
tions is reasonable, and comparable to the slightly superior forecast skill obtained 
for boreal winter conditions (Reale et al. 2008). 

3. Analysis and Simulation 

3a. Validation 

In this work we focus on two events of the SOP-3: one non-developing and one 
developing AEW, appearing quite similar in terms of the intensity of the low-level 
circulation, vertical shear and other dynamical forcings, and we investigate their 
different evolution and the different properties of the corresponding SAL intru- 
sions. Figure 1 shows a Hovmpller diagram of 700 hPa relative vorticity and rela- 
tive humidity obtained from the GEOS-5 analyses, to emphasize AEWs observed 
during the SOP-3, covering the second half of August. From Fig. 1, left panel, 
the strongest wave of the period appears on 23 August at about 5 °W, undergoing 
transition from continental to oceanic environment on the following day (hereafter 
W 1). An evident sharp strip of dry air shows the same propagation speed and am- 
plitude of the wave and is clearly associated with a Saharan Air outbreak which 
shortly precedes the wave. There is a hint of another SAL outbreak appearing on 
Fig. 1 at about 20° W and lagging of about 1 day, but this is not fully captured here 
because the outbreak occurs at higher latitude. 
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In Figure 2 the same Hovm0ller diagram is extracted from operational NCEP 
analyses (interpolated on pressure levels at 1° x 1°), for validation purposes. There 
is a very good correspondence in the vorticity fields, except that the GEOS-5 
analyses show a higher amount of detail, due to the higher resolution. As for the 
moisture, it is intrinsically a very noisy field, but the correspondence between the 
basic features is very good. 

To further confirm the legitimacy of our analyses, we also use NAMMA data 
for validation. In Figure 3, the profiles of the zonal and meridional wind com- 
ponents, obtained from the vertical grand-based soundings from the Cape Verde 
Islands, at 14.926111V, 23.4946111V and an altitude of 85m, are reported (ra- 
diosonde data are available online at: http://namma.msfc.nasa.gov). The zonal 
wind time-series clearly shows the African Easterly Jet (AEJ) located between 
700hPa and 600hPa with its amplited being modulated at time scale of two to 
seven days. The weakening of the AEJ is associated with the passage of waves, 
shown as sign inversions in the meridional wind time-series. The northerly to 
southerly change occurred at Aug 26 coincide with the passage of the positive 
vorticity over Cape Verde Islands shown in Figures 1 and 2. 

In Figure 4, zonal and meridional wind profiles extracted at the same coor- 
dinates from the GEOS-5 analyses are provided. Despite that the two data sets 
are not strictly comparable, (the GEOS-5 sequence of profiles is computed con- 
tinuously every 6 hours, whereas the plot in Figure 3 is an interpolation of daily 
profiles), the GEOS-5 analyses match very well with the observed data: the four 
episodes of maximum zonal wind on 19-21 Aug, 22Aug, 1-3 September, 5-7 
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September, and 11-13 September correspond perfectly to the observations. The 
height of the maximum wind is at about 600hPa, exactly as in the observed profiles 
and, most important, even the maximum values of zonal wind are comparable, be- 
ing on the order of 20m s _1 . 

Also a comparison of the meridional component reveal that the major wave 
passages are captured: the positive values of v, corresponding to northerly com- 
ponent (23 August, 27 August, 4 September, and 13 September) on the order of 
5 — 15m s -1 , indicating a wave passage, and present in Figure 3, are clearly re- 
produced in the GEOS-5 high resolution analyses. 

3b. Nondeveloping wave 

With the strong confidence provided by the good quality of the GEOS-5 analy- 
ses, we investigate the structure and kinematics of the major nondeveloping wave 
observed in the first part of the campaign. 

In Figure 5, 700 hPa specific humidity and 900 hPa flow are shown together 
to emphasize the relationship between two different levels in the analysis (OOz 26 
August) and across 3 forecast times (24, 48 and 72 hour forecast, corresponding 
to verification times of OOz 27, 28 and 29 August respectively). The 700hPa and 
900hPa levels correspond approximately to the vertical center and the base of the 
SAL, with the latter emphasizing the low-level circulation. Despite the strength 
of the wave, and the short-lasting apparent formation of a vertically aligned circu- 
lation evident in the analyses from lower levels up to almost 500hPa (not shown), 
W1 becomes a nondeveloping wave. Figure 5 shows that dry air at 700hPa, 
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present at OOZ 26 Aug to the west of the 900 hPa circulation center (at about 
17° IV; 21 °W), is being advected in 24 hours on the top of it (about 30 °W). After 
that, the rotating system becomes elongated and then rapidly evolves into an open 
wave (48 and 72 forecasts). This is even more evident while analyzing intermedi- 
ate time-steps (not shown). 

In Figure 6, a zonal vertical cross-section of specific humidity at 20° A" is ex- 
tracted from the GEOS-5 24-hour forecast (displayed in Fig 5), intersecting the 
SAL outbreak between 40 °W and 35°W. The signature of the SAL as a sharply 
defined ‘corridor’ of extremely dry air centered at about 38 °W (with values less 
2gKg~ 1 , surrounded by values larger than 8g/Vg _1 ) can be seen down to 800 hPa. 
Remarkable moisture gradients are present on both sides. A zonal temperature 
anomaly is obtained by subtracting the zonal mean between 80 °W and 0°: a well- 
defined thermal dipole, stronger than any other anomaly in the range of longitudes 
selected, (represented by a warm anomaly between 800hPa and 400hPa and a cool 
anomaly between 825 hPa and the surface) can be seen in correspondence to the 
dry tongue. Since the cross-section cuts across the SAL intrusion, it appears that 
temperature, in the core of the SAL, is up to 3°C warmer than the surroundings at 
the same latitude. Most important, in the section there is also evidence of a thin 
dry air intrusion being entangled in the circulation at about 48° — 51 °W . This 
intrusion is not associated to a thermal dipole. A step-by-step streamline analysis 
compared with microwave geostationary imagery (not shown) show quite clearly 
that this is a fine intrusion of midlatitude air, produced by a large anticyclone in 
the Atlantic, dry because of subsidence. The GEOS-5 0.25° forecasts keep these 
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two air masses well separate. 

In Figure 7, for verification purposes, the same figure is extracted from the full- 
resolution NCEP operational analyses in model sigma levels (0.5 corresponding 
roughly to 500 hPa). Despite small-scale differences, the NCEP analysis confirms 
the thermal structure depicted in the GEOS-5 24-hour forecast and the presence 
of a very well-defined dry intrusion at about 35° — 40 °W. Most remarkable is 
the presence of a similar thermal dipole seen in Figure 3. A positive anomaly 
of up to about 3°C, even stronger than the GEOS-5 forecast, is present in the 
mid-troposphere, and a corresponding cool anomaly in the low moist layer. Both 
anomalies are obtained, as in Figure 3, by subtracting the 80°W — 0° mean. This 
is in partial agreement with Dunion and Velden, (2004). The relatively new aspect 
of this analysis is that also a negative value appears in the moist low-level layer at 
the base of the column. Finally, it is important to notice that the thin midlatitude 
dry air intrusion, clearly seen in the GEOS-5 analyses and maintained in the 24- 
hour GEOS-5 forecasts as a separate air mass (Figure 6), is not represented in the 
NCEP analyses (Figure 7), possibly being merged with the surrounding air mass. 

3c. Developing wave 

The Hovmpller computed in Figure 1 shows a strong wave moving from about 
10 °W to 40°!F between OOz 10 and 15 September (hereafter W2). Saharan Air 
precedes the wave: however, in contrast with Wl, values of low humidity are less 
extended, suggesting that the SAL intrusion was not at strong. This is the wave 
that is associated with the system eventually named as Helene. 
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In Figure 8, the GEOS-5 forecast for the genesis of Helene is shown. The 
900 hPa circulation shows a clearly defined vortex (at about 25°W and 12 ° N 
in the analyses) progressing westward and then recurving northwestward, being 
entangled in corresponding high levels of 700 hPa moisture. The displacement 
error of the predicted center with respect to the observed position is smaller than 
2° at all times (not shown). 

In Figure 9 a zonal cross-section at 13°iV (similarly to what was done to in- 
vestigate W1 in Figure 6) is produced, intersecting both the center of the system 
W2 (which is a precursor of Hurricane Helene), and the SAL intrusion. The ver- 
tically aligned vorticity column at about 30° — 33 °W, stretching from the surface 
to the the tropopause, is the signature of the Tropical Storm, named Helene at 
OOz 14 September 2006 (Brown, 2006). A positive temperature anomaly of more 
than 2° at about 35 — 45 °W in the lower midtroposphere is associated with the 
same Saharan air outbreak which in the Hovmpller in Figure 1 reaches 40 °W on 
14 September. The anomalies seen within the vorticity column associated to He- 
lene are not associated to the SAL but are due to the tropical storm development 
(a hint of a warm core above 500hPa and some cooling possibly due to rainfall 
evaporation in the lower levels). 

Comparing this cross-section with the one relative to the nondeveloping W 1 
(Ligure 6), two prominent differences can be seen: there is minimal cool anomaly 
in the lowest levels, and the dry air appears more diluted with less sharp horizontal 
gradients. 

In Ligure 10, for validation purposes, we plot the same vertical zonal cross- 
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section seen in Figure 9, but this time from the full-resolution NCEP analyses in 
sigma levels. The vertical column of vorticity present at about 30° — 3 3° IF is 
the analyzed signature of Tropical Storm Helene, and a warm anomaly of up to 
3 °C at about 40°, corresponding to the SAL intrusion, is also present in the mid- 
lower troposphere. Despite the much weaker NCEP representation of Helene with 
respect to the GEOS-5 (the latter being closer to observations since Helene at that 
time was already a named system) the NCEP and GEOS-5 systems agree well on 
the overall scenario of a SAL intrusion preceding Helene, but not as strong as the 
case seen in W1 . 

4. Discussion 

4a. Possible role of the SAL 

W1 and W2 do present some similarities: among them, vertical shear, computed 
as difference between the vertically integrated 100-200hPa and 800-900hPa wind 
speed, does not reach 10m s -1 . Other environmental conditions, such as sea sur- 
face temperature and low-level vorticity appear very favorable in both W1 and 
W2 case (not shown); however only the latter underwent development becoming 
Helene. The only difference appears to be the intensity of the Saharan Air intru- 
sion depicted in analyses and GEOS-5 forecasts. In the model, the temperature 
dipole associated with the SAL affecting W 1 can be followed at each time step 
and can be considered a possible cause of suppression (increasing static stability 
in the wave) together with the SAL dryness. In the precursor of Helene, while dry 
air being present especially on its northern and western quadrants, the low-level 
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negative anomaly is minimal or absent, and therefore the effect on static stability 
is smaller. 

The positive anomaly can be simply attributed to the signature of warm air 
originated over the Sahara but the cool anomaly in the lower levels does not have 
any plausible explanation relying on transport only. There is no source of lo- 
calized cooler temperatures at that latitude, away from landmass and in a very 
homogeneous marine tropical environment. At this time a possible explanation is 
that the low-level cool temperatures are at least partly an indirect evidence of dust 
amount. The thermal effect of Saharan mineral dust is a net reduction of down- 
wave shortwave radiation in the near- surface levels, and a heating in the lower 
midtroposphere, corresponding to the core of the SAL. 

To support this possibility, we present in Fig. 11 the aerosol optical depth 
measured from the Moderate Resolution Imaging Spectroradiometer (MODIS) is 
compared between the periods of 26-29 August 2006 and 13-16 September 2006. 
These periods encompass the snapshots of the circulations in Figures 5 and 8. 
Notwithstanding some contamination of MODIS data by clouds, the air advected 
by the large-scale flow in W1 circulation, especially on its northern side, is likely 
to contain a larger amount of dust. The different amount of dust must have a dif- 
ferent impact on the thermal structure of the atmosphere. It appears that the high- 
resolution analyses and forecasts can represent these different thermal structures 
and that the GEOS-5 model initialized by such analyses can retain it for 24-72 
hours advecting it into the circulation and producing a realistic cyclone dissipa- 
tion. In the case of the GEOS-5 forecast, the nature of the finite-volume dynamics 
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(Lin, 2004) is such that is a particularly suitable tool to generally maintain sharp 
gradients by minimizing unrealistic diffusion processes. The finite-volume dy- 
namics has been shown to be very efficient in the mid latitudes where localized 
temperature gradients associated with sharp fronts can be very realistically simu- 
lated and maintained. This work documents that the same skill can be very useful 
also in the tropics when dealing with Saharan Air. 

4b. Critical role of data coverage and model resolution 

In this work, as stated in the introduction, two sets of forecasts were produced 
from the GEOS-5 high quality analyses, at resolution of half and a quarter of a 
degree. The most prominent result is that the forecasts issued at a half degree, 
when compared to the quarter degree, are not as suitable to retain all the infor- 
mation present in the analyses, and, especially with respect to the SAL, tend to 
irreparably blur the sharpness of its boundaries. Also the thermal information is 
sometimes diluted and spread unrealistically. To provide evidence of this prob- 
lem, in Figure 1 1 shows the same section of Figure 6 (which was obtained from 
the 24-hour forecast at 0.25°) but from the corresponding forecast at 0.5°. The 
most prominent difference is that there the SAF filament is much broader and 
with less sharply-defined borders. Moreover, there is no clear separation between 
the midlatitude dry air at about 50 °W and the SAF (at about 35° — 42°) as seen 
in Figure 6. In addition, the cool anomaly in the lower level is substantially less 
confined and spread in the lower resolution forecasts. Finally, Figure 6 shows a 
moist column at 25° — 30° W which is associated with W 1 ; this feature is basically 
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absent from the half degree forecast. 

It is important to stress that it is not only the resolution of the model or of the 
data assimilation system which can suppress the fine structure of the SAL and its 
thermal properties. On data poor area, such as the Sahel and the eastern tropical 
Atlantic the optimal use of satellite information, far for being implemented in any 
operational system, is probably as important as the resolution. 

One of the purposes of the GEOS-5 development was to contribute to the 
creation of a new set of analyses, the Modem Era Retrospective Analysis for Re- 
search and Applications (MERRA), documented by Bosilovich (2008) by making 
a much more extensive use of all satellite information. The GEOS-5 used in this 
work is almost the same version adopted for MERRA. The difference between the 
NCEP operational analyses and the GEOS-5 analyses (which are of comparable 
resolution) evident for example between Figures 9 and 10 in terms of sharpness of 
moisture boundaries, and the intensity of Helene, are probably due to the different 
assimilation of satellite data. 

5. Conclusion 

This work investigates AEW development during the the SOP-3 phase of the 
NAMMA campaign with the aid of the NASA GEOS-5 data assimilation and 
forecasting system. A data assimilation run is performed to produce a high qual- 
ity set of global analyses which make heavy use of satellite information. From 
these improved analyses, two sets of 30 5-day forecasts at two different resolution 
are produced. Forecasts are validated against operational NCEP analyses. 
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A very strong non-developing wave is selected and is compared with the wave 
that eventually developed into Hurricane Helene. Analyses and forecasts docu- 
ment the presence of a strong temperature dipole associated with the Saharan air 
intrusion. This dipole is advected into the circulation of the wave, and may have 
contributed to inhibit further development by increasing static stability. No such 
dipole is found for the Saharan air that intruded in Helene’s precursor. MODIS 
optical thickness show that the SAL outbreak of the first case contained a higher 
amount of dust. The lower tropospheric cooling associated with the strong Saha- 
ran air outbreak suggests that the high resolution global simulations and analyses 
can capture not only the SAL intrusion but also possibly part of the thermal effect 
consequent to downward shortwave reduction caused by large amounts of Saharan 
dust. When this thermal structure is well-represented in the initial conditions, a 
high-resolution global model can have more skill in predicting the impact of SAL 
on a tropical cyclone. 

At this time, the role of dust on tropical cyclogenesis is still subject to a vigor- 
ous debate and contrasting evidence have been presented by Dunion and Velden 
(2004) in which a cyclogenesis suppression is suggested, and Jenkins et al. (2008) 
which provide evidence of rain band invigoration attributed to dust. 

It is clear that the details of possible radiative effects by dust on cyclogene- 
sis will have to be investigated by a high-resolution GCM with interactive dust 
aerosol. This capability is not yet available anywhere within the scientific com- 
munity, but this work aims to provide a set of important needed requirements to 
even plan such research. 
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In particular, we have emphasized how the resolution is a critical parameter 
in being able to represent with some degree of confidence the fine sharpness of 
SAL intrusions in tropical circulations. No simulation can be reasonably hope to 
investigate the complexity of the mechanisms involved with interaction of SAL 
and tropical cyclones at a resolution lower than a quarter of a degree. This is a 
demanding standard for global modeling but the evidence provided in this paper 
suggests that it is an inescapable necessity. 

Second, in this work it is shown that even at the respectable resolution of the 
current operational analyses, inadequate use of satellite data, especially on data 
poor areas, may be a serious limitation. On the other hand, the density of satellite 
information adopted and, most important, the quality of the coverage, may bring 
substantial improvement to our understanding of processes associated to tropical 
cyclogenesis. 

Acknowledgments Authors acknowledge support from Dr. Ramesh Kakar, NASA 
Headquarters, through the NAMMA Project, and use of NASA High-End Com- 
puting resources. 
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Figure 1: Hovm0ller of 700 hPa relative vorticity (s -1 , left panel) and relative humidity 
(right panel) from the GEOS-5 operational analyses, latitudinally averaged (12° — 20°iV), 
covering the SOP-3 period from 15 August to 15 September. 
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Figure 2: Hovm0ller of 700 hPa relative vorticity (s _1 , left panel) and relative humidity 
(right panel) from the NCEP operational analyses, latitudinally averaged (12° — 20° iV), 
covering the SOP-3 period from 15 August to 15 September. 
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Figure 3: Sequence of the observed profiles of the zonal (above) and meridional (be- 
low) components of the wind, taken daily at 12Z at the Cape Verde station (Latitude 
14.926111°iV, Longitudie 23.494611 °W, and altitude 85m) during the SOP-3 campaign. 
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Figure 4: Same as Figure 3, but extracted from the GEOS-5 analyses: sequence of the 
observed profiles (at each synoptic time) of the zonal (above) and meridional (below) 
components of the wind, taken at the same latitude and longitude of the sation. 
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Figure 5: GEOS-5 (at 0.25°): 700 hPa specific humidity (g Kg -1 ) and 900 hPa wind 
(streamlines) in the GEOS-5 analyses (upper left) for OOz 26 August, and relative to the 
24, 48 and 72 hour forecasts for OOz 26, 27 and 28 August. 


30 



100 


GEOS5 Lat20N IN00z26Aug06 24hFc VER27Aug 



Figure 6: GEOS-5 (at 0.25°): zonal vertical cross-section of specific humidity {g Kg" 1 , 
shaded) and temperature anomaly (°C, red/blue contour, positive/negative, subtracting 
the zonal mean between 80 °W and 0°) at 20°A r for OOz 27 August, 24 hour forecast 
initialized at OOz 26 August. W 1 located at about 25 — 30 °W. 
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Figure 7 : As Figure 6, but extracted from NCEP full-resolution analyses in model levels 
for validation purpose: zonal vertical cross-section of specific humidity ( gKg~ l , shaded) 
and temperature anomaly (C, red/blue contour, positive/negative, subtracting the zonal 
mean between 80°FF and 0°) at 20 °N for OOz 27 August. W 1 located at about 25 — 30 °W 
The vertical dimension is only approximately comparable with Figure 3 since the spacing 
between model levels and pressure levels is different. 
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Figure 8: GEOS-5 (at 0.25°): 700 hPa specific humidity ( gKg _1 , shaded) and 900 hPa 
wind (streamlines) in the Analyses (upper left) for OOz 13 September, and relative to the 
24, 48 and 72 hour forecasts for OOz 14, 15 and 16 September. 
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Figure 9: GEOS-5 (at 0.25°): zonal vertical cross-section of specific humidity {g Kg~ l , 
shaded), temperature anomaly ( °C , red/blue contour, positive/negative, subtracting the 
zonal mean between 80 °W and 0°) and cyclonic vorticity (darker red, contours at 1 , 2, 3 x 
10 _4 s~ 1 at 20°iV for OOz 14 September, 24 hour forecast initialized at OOz 13 September. 
The developing TS Helene, identifiable by the vertical column of strong cyclonic vorticity, 
is located at about 2>2°W . 
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Figure 10: As Figure 9, but extracted from NCEP full-resolution analyses in model levels 
for validation purpose: zonal vertical cross-section of specific humidity ( g Kg ~ 1 , shaded), 
temperature anomaly (°C, red/blue contour, positive/negative, subtracting the zonal mean 
between 80 °W and 0°) and cyclonic vorticity (darker red, contours at 1, 2, 3 x 10 -4 s -1 
at 20 °N for OOz 14 September. The developing TS Helene, identifiable by the vertical 
column of strong cyclonic vorticity, is located at about 32 °W. 
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Figure 1 1 : Four-day mean MODIS aerosol optical depth for (a) 26-29 August 2006 and 
(b) 13-16 September 2006. The MODIS data show much larger amount of dust for the 
first wave, which is the non-developing one. 
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Figure 12: As Figure 6, but produced from GEOS-5 at the lower resolution of 0.5° zonal 
vertical cross-section of specific humidity ( g Kg 1 , shaded) and temperature anomaly 
(°C, red/blue contour, positive/negative, subtracting the zonal mean between 80° W and 
0°) at 20 ° N for OOz 27 August, 24 hour forecast initialized at OOz 26 August. 
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